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Abstract. We analyze ISO archive data of the nearby bright emission-line star 51 Oph, previously classified as a proto- 
planetary system similar to j3 Pic. The infrared spectrum reveals the presence of gas-phase emission bands of hot (~ 850 K) 
CO, CO2, H2O and NO. In addition to this, partially crystalline silicate dust is present. The solid-state bands and the energy 
distribution are indicative of dust that has formed recently, rather than of debris dust. The presence of hot molecular gas and the 
composition of the circumstellar dust are highly unusual for a young star and are reminiscent of what is found around evolved 
(AGB) stars, although we exclude the possibility of 51 Oph belonging to this group. We suggest several explanations for the 
nature of 51 Oph, including a recent episode of mass loss from a Be star, and the recent destruction of a planet-sized body 
around a young star. 

Key words. Circumstellar matter - Stars: emission-line - Stars: evolution - Stars: 51 Oph - Infrared: Stars 



1. Introduction 

51 Ophiuchi (HR 6519) is a bright (V = 4.78), nearby (d = 
131 pc) emission-line star. Although the star has traditionally 
been classified as B9.5Ve (Buscombe 1963; Slettebak 1982; 
Jaschek & Jaschek 1992; Dunkin et al. 1997), recent work by 
Gray & Corbally (1998) show 51 Oph to be of spectral type 
AOII-IIIe, more in line with its position to the right of the 
main-sequence in the Hertzsprung-Russell diagram (van den 
Ancker et al. 1998). For the current generation of ground-based 
telescopes, 5 1 Oph is a point-source at visual to mid-infrared 
wavelengths (Merkle et al. 1990; Lagage & Pantin 1994). 

51 Oph exhibits a normal UV and visual energy distribu- 
tion, but IRAS and subsequent ground-based follow-up obser- 
vations revealed an unusually large infrared excess longward 
of 2 ^m, interpreted as being due to hot (up to 1000 K) cir- 
cumstellar dust (Cote & Waters 1987; Waters et al. 1988). The 
peculiarity of the energy distribution of 5 1 Oph when compared 
to those of Be, Ae and A-shell stars led Waters et al. (1988) to 
suggest that it may be a candidate proto-planetary disk system, 
similar to f3 Pictoris. This suggestion was further strengthened 
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by Grady & Silvis (1993), who found evidence for the presence 
of variable columns of accreting gas, similar to those found 
around (3 Pic. 

The presence of dust in the 5 1 Oph system was unambigu- 
ously established in 1993, when Fajardo-Acosta et al. detected 
a prominent 10 //m silicate feature in emission. Their results 
were subsequently confirmed from the ground by Lynch et al. 
(1994), Walker & Butner (1995) and Sylvester et al. (1996) and 
by Waelkens et al. (1996) usingthe Infrared Space Observatory 
(ISO; Kessler et al. 1996). Lecaveher des Etangs et al. (1997) 
found cold neutral C (iV = 5 x lO^^ cm^^, T = 20 K) in 5 1 Oph, 
yielding a C l/dust ratio similar to that derived for f3 Pic. Since 
C I has a very short lifetime, it must be continuously replen- 
ished, providing evidence for the existence of evaporating bod- 
ies in 51 Oph, similar to those infeiTed around f3 Pic. 

However, the analogy with the /3 Pic system is not com- 
plete. 5 1 Oph possesses strong Ha emission, absent in /? Pic. 
The shape of the spectral energy distribution suggests that only 
very few large dust grains are present, unusual for a protoplane- 
tary disk, but in agreement with what is found in Vega-type sys- 
tems. However, the magnitude and temperature of the infrared 
emission is more reminiscent of disks surrounding younger 
stars. These indicators, together with its position to the right 
of the main-sequence, lead some authors to consider 5 1 Oph to 
be related to the Herbig Ae/Be stars (e.g. van den Ancker et al. 
1998;Meeus etal. 2001). 
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Fig. 1. ISO-SWS spectra of 51 Oph. The main figure shows the AOT 01 spectrum, whereas the insets show the AOT 06 data. 
The grey Unes show a Kurucz (1991) model for the AO primary (dark) and the spectrum (light) obtained by subtracting this 
photosphere model from the observed spectrum. 



In this letter we will re-analyze archive ISO spectroscopy 
of 5 1 Oph. We will show that its 4-8 /xm spectrum, not shown 
by Waelkens et al. (1996), reveals the presence of hot circum- 
stellar molecular gas. This is highly unexpected for a proto- 
planetary or debris disk system. We investigate several different 
explanations for the remarkable composition of the circumstcl- 
lar material in 5 1 Oph and conmient on those. 

2. Observational data 

ISO observed 5 1 Oph with the Short Wavelength Spectrometer 
(SWS; de Graauw et al. 1996) in mode "AOT 01" (full spec- 
tral scan) and twice in mode "AOT06" (deeper scans over a 
limited wavelength range). In addition, a Long Wavelength 
Spectrometer (LWS; Clegg et al. 1996) full grating scan ("AOT 
LOl") was taken. The data were reduced in a standard fashion 
using calibration files corresponding to ISO off-line processing 
software (OLP) version 7.02, after which they were corrected 
for remaining fringing and glitches. The resulting SWS spectra 
are shown in Fig. 1 . 

The spectra observed by ISO are in good agreement with 
the ground-based iV-band spectroscopy of 51 Oph by Fajardo- 
Acosta et al. (1993), Walker & Butner (1995) and Sylvester 
et al. (1996). However, the more complete wavelength cov- 
erage, and the higher spectral resolution and signal to noise 
reveal a wealth of structure hidden in the ground-based data. 
Shortward of 4.1 /xm, the spectrum consists of a smooth con- 
tinuum, dominated by the photosphere of the central star. A 
broad absorption due to Br/3 can be seen around 2.63 /xm. 



The Bra line at 4.05 ^m is filled in with emission (line flux 
2.5 X 10~^^ W m~^), confirming the emission-line classifica- 
tion of 51 Oph. 

A number of strong spectral features due to gas-phase 
molecules are present in the 4.1-7.5 /xm spectral range. We 
recognise the 4.2-5.4 fim fundamental vibritional band (v = 
1-0) of CO, the P and R branches of the 1^2 bending mode of 
H2O in the 5.2-7.5 i^m range, the 4.1-4.4 /xm 1/3 band of CO2 
and the fundamental vibritional band of NO (At; = 1) around 
5.1-5.9 fim. To further quantify the gas-phase emission lines, 
we have compared the continuum- subtracted 5 1 Oph spectrum 
to a model for gas-phase emission using a single layer, plane 
parallel LTE code (Cami et al. 2000). A satisfactory fit to our 
data can be obtained with T^oi. = 750-900 K, iV(CO) = lO^"- 
10^2 cm-2. The best fitting model has T^oi. = 850 K, 7V(H20) 



cm--', NC^C02) = 3 X 10^*' cm^^ 7V(^^CO) = 



= 4 X IQis 

3 X lO^i cm-^, NC^CO) = 4 X 10^^ cm^^ and N^^NO) 
1 X 10^^ cm~^. It is shown in Fig. 2. We stress that this fit is 
not unique and that several other combinations of temperature 
and column density also give a satisfactory fit. However, we 
can exclude a large abundance of ^^CO. The presence of NO 
is surprising but is needed in order to fit the emission structure 
near 5-5.5 iJ.m. 

Longward of 8 /xm, the spectrum is dominated by solid- 
state emission from 0-rich dust. A strong amorphous silicate 
feature is present, peaking at 10.3 /ttm. As was already noted 
by Fajardo-Acosta et al. (1993) and Sylvester et al. (1996), the 
long-wavelength shoulder of this feature consists of an almost 
hnear part extending up to 15 /zm. We note that this behaviour 



M.E. van den Ancker et al.: The composition of circumstellar gas and dust in 51 Oph 



3 



5n • • \ • • • • \ • • • • \ • • • • \ • ^ 



- Total fit ! 
- Observation : 




4.5 5.0 5.5 6-0 

A [/im] 



Fig. 2. Continuum subtracted 4-6 fjia spectram of 51 Oph 

(heavy soHd Hne) compared to a model for gas-phase emis- 
sion of H2O, CO, CO2 and NO at 850 K with iV(H20) = 
4 X IQi* cm-2, iV(i2C02) = 3 x 10^^ cm'^, iV(i2cO) = 
3 X 10^1 cm-2, 7V(i3CO) = 4 X 10^^ cm-2, and 7V(14nO) 
1 X 10^^ cm~2 (grey line). Contributions of the different gas- 
phase components are indicated by the shaded area. 

is incompatible with a purely silicate origin of this feature. A 
secondary feature, peaking around 11. 5 /xm and with a width of 
~ 3 /xm must be present as well to be able to explain the shape 
of the 8-15 fim emission complex. A very broad 10 /im band is 
also seen in the enigmatic object 77 Car (Morris et al., in prepa- 
ration). 

A second broad bump around 18 /im is present and may 
be due to amorphous silicates. In the SWS AOT 06 spectrum 
a narrower (FWHM ^ 0.4 /itm) peak at 19.2 /xm, similar to 
that seen in stars with emission due to crystalline silicates (e.g. 
Molster et al. 1999), is superimposed on the broad amorphous 
sihcate feature. 51 Oph was not detected in the LWS scan. 

3. Discussion and Conclusions 

The presence of spectral features in 51 Oph due to hot gas- 
phase molecules is highly surprising. None of the 47 Herbig 
Ae/Be stars for which material is present in the ISO data 
archive shows a 4-8 fim emission complex resembling the one 
found in 51 Oph. However, such a complex due to emission 
from hot molecular gas is commonly found in the ISO spectra 
of O-rich evolved stars on or near the Asymptotic Giant Branch 
(AGE; e.g. Yamamura et al. 1999, Sylvester et al. 1999). 

A further inspection of the ISO data archives confirms our 
suspicion that the spectrum of 5 1 Oph is highly unusual for a 
protoplanetary system: none of the young stars for which spec- 
troscopic data is present shows the 11.5 /xm shoulder to the 
amorphous silicate feature found in 5 1 Oph, including the iso- 
lated Herbig Ae/Be stars, of which 5 1 Oph is beheved to be a 
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Fig. 3. Continuum subtracted 7-16 /xm spectrum of 51 Oph 
(black sohd line) compared to continuum- subtracted ISO SWS 
spectra of the Herbig Ae star AB Aur (van den Ancker et al. 
2000; dark grey curve) and that of the symbiotic star V835 Cen 
(Mikolajewska et al. 1997; light grey line). 

member (Meeus et al. 2001; Bouwman et al. 2001). However, 
comparison of the observed 10 /xm profile with archive ISO- 
SWS data of the symbiotic star V835 Cen (Fig. 3) shows that 
the 10-15 /im complex in 51 Oph is similar to the one found 
in this highly evolved system. In fact, the presence of a long- 
wavelength shoulder to the 10 /im feature is a common feature 
in the spectra of O-rich AGB objects (e.g. Tielens et al. 1998). 
It is not observed in other types of stars. 

The question arises what is the origin of the circumstellar 
gas and dust in 5 1 Oph. In view of the peculiar molecular spec- 
trum (indicative of a large column of hot gas) and the deviating 
10 /xm silicate band, as well as the sharply dropping mid-IR to 
millimeter continuum (pointing to a lack of large, cold grains 
typical for proto-planetary disks) it is interesting to consider 
the possibiUty that the circumstellar dust in 51 Oph is "fresh", 
i.e. recently formed rather than accreted from the interstellar 
medium and processed in a proto-planetary disk. 

One possible explanation for the presence of hot gas and 
newly formed dust in 5 1 Oph would be to assume that 5 1 Oph 
is a highly evolved rather than a relatively young system. This 
assumption agrees well with the above main-sequence position 
of 51 Oph in the HRD. The optical spectrum of 51 Oph, as 
well as its luminosity, show that a single post-main sequence 
star can not have created the large amount of circumstellar ma- 
terial we observe. However, the presence of a cool post-main 
sequence companion may explain the observed properties of 
the system. We note that a companion to 51 Oph may have 
been detected: Buscombe (1963) found the radial velocity of 
51 Oph to be variable, which could be due to the reflex motion 
of the companion. 

By comparing the optically visible AO star's TcH of 
10,000 K and L/Lq of 260t^g Lg with the post-main se- 
quence evolutionary tracks by Schaller et al. (1992), we derive 
a mass of 3.8 Mq and an age of 2 x 10^ years for the primary. 
Using the ISO data and the Hipparcos distance of 1311^3 pc 
towards 51 Oph we compute an infrared excess of 7.1'^2a L© 
for 5 1 Oph. Assuming that both stars in our putative binary sys- 
tem formed at the same time, a lower mass companion should 
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stil be on the main-sequence, limiting its mass to < 2 Mq if its 
luminosity is not to exceed the observed infrared excess. This 
lower mass for the more evolved star in the system is at odds 
with the mass derived for 5 1 Oph itself, unless significant in- 
teraction between both stars occurred. In addition, it is highly 
unUkely that an evolved star of less than 7 Lq shows such a 
large amount of circumstellar matter. 

Another possibility would be to assume that the companion 
is more massive than the optically visible AO star. Any scenario 
in which the companion would be a red giant or AGB star of 
comparable mass to the primary is impossible, since such a star 
would have a luminosity that greatly exceeds 7 L0 . However, 
if the companion would be sufficiently massive 7 Mq), it 
would already have evolved to become a massive white dwarf 
of very low luminosity. In this scenario, the material observed 
around 5 1 Oph today could be caught in remnants of material 
from the AGB wind of its companion. However, the lack of a 
significant amount of ^^CO in the composition of the circum- 
stellar material argues against this scenario. We conclude that 
the presence of an evolved companion in 51 Oph is unlikely. 

The rapidly rotating (vs'mi = 270 km s~^; Dunkin et al. 
1997), emission-Une character of 51 Oph has led to an asso- 
ciation with the Be stars, who share many of the properties of 
51 Oph. However, as pointed out by Waters et al. (1988), cir- 
cumstellar dust is highly unusual in Be stars. The presence of 
double-peaked Ha emission in 51 Oph does indicate the pres- 
ence of a gaseous disk, and the width of the line suggests that 
this gas is close to the star. It is reasonable to assume that the 
dust is also located in this disk, at a larger distance from the star. 
In Be stars, the disk can either be the result of (time- variable) 
mass loss, or from mass accretion from an evolved companion. 
However, as pointed out above, the presence of an evolved cool 
companion is difficult to understand given the luminosity of the 
dust/gas emission in the infrared. If one accepts that the circum- 
stellar material is due to a recent period of high mass loss from 
the Be star into a disk, the presence of dust suggests that the 
amount of gas ejected must have been high enough to shield 
the UV radiation field of the central star to allow for dust for- 
mation. We have fitted the ISO-SWS spectrum using the dust 
radiative transfer code MODUST (Bouwman & de Koter, in 
preparation) and find a dust mass of 4 x 10~^ Mq. This is a 
factor of 5 higher than the (distance-corrected) value given by 
Fajardo-Acosta et al. (1993). In this scenario the late-main se- 
quence nature of 51 Oph may be related to the ejection of a 
large amount of inass, since rapidly rotating early type stars are 
believed to move closer to their breakup velocity as they evolve 
(e.g. Langer et al. 1999; Langer 2000). 

A more exotic possibility to explain the properties of 
5 1 Oph would be to infer that we are seeing the aftermath of 
a recent event such as the collision of two gas-rich planets or 
the accretion of a solid body as the star increases its size at the 
end of its main-sequence hfe. New dust could form at the site 
of the evaporation of the solid body, explaining both the high 
column of hot gas, the apparent small dust particle sizes and 
the composition of the silicate dust. 

As was already pointed out by previous authors (e.g. Herbig 
1994) the problem of separating young stars from their evolved 
counterparts is non trivial. The case of 5 1 Oph presented in this 



letter may serve once more to illustrate this difficulty. Even 
with the extensive data present on this enigmatic object, its 
evolutionary status remains unclear. New data that is able to 
resolve the system spatially, as well as a more thorough in- 
vestigation of the suspected radial velocity variations, may be 
needed to shed more light on its nature. 
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